Biaxial Deformations of Rubber: Entanglements or Elastic Fluctuations? 
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The classical theory of rubber elasticity fails in the regime of large deformation. The underlying 
physical mechanism has been under debate for long time. In this work, we test the recently pro- 
posed mechanism of thermal elastic fluctuations by Xing, Goldbart and Radzihovsky [l( against the 
biaxial stress- strain data of three distinct polymer networks with very different network structures, 
synthesized by Urayama 2] and Kawabata [3j] respectively. We find that both the two parameters 
version and the one-parameter version of the XGR theory provide satisfactory description of the 
elasticity in whole deformation range. For comparison, we also fit the same sets of data using the 
slip-link model by Edwards and Vilgis with four parameters. The fitting qualities of two theories 
are found to be comparable. 

PACS numbers: 62.20.de, 61.41. +e 



I. INTRODUCTION 

The classical theory of rubber elasticity Q properly 
captures the essence of entropic elasticity of polymer net- 
works. As it is well known, however, it fails to predict 
accurately the stress-strain relation at large deformation 
regime. Understanding of the mechanism of this fail- 
ure and attempting to improve it has remained a central 
problem in field of polymer science in the last half cen- 
tury. 

The classical theory is based on two basic assump- 
tions: 1) The probability distribution of the end-to-end 
displacement of polymer chains are Gaussian; 2) the to- 
tal entropy of the network is the sum of entropy of each 
chain. Regarding the assumption 1), while it is well 
known that polymer chains are Gaussian, i.e. ideal, in 
a polymer melt, wether this remains true after the sys- 
tem being crosslinked is not so clear. It is however the 
fundamental assumption of almost all entanglement the- 
ories that this assumption is incorrect due to topological 
constraints imposed by entanglements. Exact treatment 
of entanglement effects is strictly speaking a complicated 
topological problem, and does not seem possible at this 
stage. Hence the possibility that entanglements do not 
invalidate the classical theory in the limit of long chain 
can not be excluded. Nevertheless, different researches 
have made distinct assumptions about the effects of en- 
tanglements in rubber elasticity, and have developed dif- 
ferent versions of entanglement theory. These include the 
Slip-Link Model by Edwards and Vilgis [5], the diffused- 
constraint model by Kloczkowski, Mark, and Erman [6|, 
the localization model by Gaylord and Douglas 0, [8[ , as 
well as the slip-tube model by Rubinstein and Panyukov 

i. 

The assumption 2) seems much less challenged, but 
can be easily seen to be invalid. It is obvious that the 
entropy of network junctions is nonzero and is not in- 
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eluded in the classical theory. James and Guth [10] has 
shown that, for phantom networks, this part of entropy 
is independent of strain deformation, and therefore has 
no effect on the elasticity. Historically, this result may 
has been seen as an evidence supporting entanglement 
theories. Real rubbery materials are, however, nearly 
incompressible, with Poisson ratio 1/2. Junction fluctu- 
ations are essentially different in phantom networks and 
in incompressible networks. This difference was made 
clear in a recent work pJ] by Xing, Goldbart, and Radzi- 
hovsky (hereafter referred as XGR). At long length scale, 
these fluctuations can be simply modeled as thermal elas- 
tic fluctuations, i.e. phonon fluctuations, of an isotropic 
solid. In incompressible solids, longitudinal phonons are 
infinitely massive and therefore do not renormalize the 
elasticity. It turns out that the fluctuations of trans- 
verse phonon modes depend sensitively on the macro- 
scopic strain deformation, and are quantitatively compa- 
rable to the effect of single chain entropy addressed by 
the classical theory. This strong interplay between long 
wave length thermal fluctuations and elastic properties 
is not peculiar to rubber at all. In fact it has been one 
of the main theme of modern condensed matter physics, 
and has been analyzed in many different systems, such 
as polymers, membranes, and liquid crystals. 

Both entanglement theories and XGR theory fit well 
to uniaxial stress-strain data. To distinguish these the- 
ories, therefore, it is important to test them aganist 
data of more general biaxial deformations. [15| Recently 
Urayama et. al. [2] tested five different versions of entan- 
glement theories using biaxial stress-strain data of a well- 
characterized end- linked networks. The slip-link model 
by Edwards and Vilgis, hereafter referred as EV, was 
found to provide the best fit. In a more recent workfll], 
Hansen, Skov, and Hassager tested XGR theory and a 
version of slip-link model against the unaxial stress-strain 
data by Mark et al^ and concluded that XGR and slip- 
link model fit the uniaxial data with comparable quality. 
In this work, we use two sets of biaxial data by Urayama, 
as well as another set of earlier data by Kawabata et. al. 
[3| to test the XGR's elastic fluctuation theory, and com- 



pare it with EV's slip-link theory. 



II. THEORIES 



of polymer chains. Ci is the polymer chain density, while 
C2 is the slip-link density. This theory reduces to the 
classical theory when both 77, a vanish. 



Let an incompressible isotropic solid be biaxial de- 
formed. A mass point x is deformed to x' = X- x, where 
the deformation gradient matrix A has three eigenvalues 
Ai, A2, A3 in three orthogonal directions. The incompress- 
ibility constraint imposes the relation: A3 = l/(AiA2); 
therefore only two of these three variable can be inde- 
pendently controlled. If we impose deformations Ai, A2 
along two axes and let the third axes relax freely, the 
two components of nominal (engineer) stress (measured 
relative to the undeformed area) are given by Q 

'^ --'^ (2.1) 
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with A3 treated as a function of Ai , A2 . The nominal 
stress (J3 along the third axis is zero, since the third di- 
mension can relax freely. 

According to the XGR theory, the total elastic free en- 
ergy density is characterized by two tunable parameters: 

f{\) = ^Trg + fi,{logq-g-^-q)^. (2.2) 

Here the average ( • )q is over the orientation of the 3d 
unit vector g, and g = A A is called the metric tensor. 
The first term is just the classical rubber elasticity theory 
and is due to the entropy of polymer chains. The param- 
eter /io has the physical interpretation of the effective 
density of polymer chains that contributing to the net- 
work elasticity. This value however may receive substan- 
tial contribution from entanglement effects, and therefore 
may be significantly different from the real chain density. 
[l6j The second term is the lowest order correction due to 
elastic fluctuations; /ii is half the total number of phonon 
modes in unit volume. Two parameters /io^/^i are there- 
fore expected to be of the same order of magnitude. In 
reality, however, their best-fitting values are always found 
to be very close. There maybe relevant physics underly- 
ing this intriguing coincidence. 

The slip-link model by Edwards and Vilgis |5l] models 
the effect of entanglement between two chains as a slip- 
link. The free energy density [17] is characterized by four 
parameters: 
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where S = A^+A2+A3, 77 is the slippage parameter, being 
zero if the entanglement behaves completely as a cross- 
link; the parameter a characterizes the finite extensibility 



III. EXPERIMENTAL SYSTEMS 

The two systems synthesized by Urayama et. al. 
[l2[ are end-linked poly(dimethylsiloxane) (PDMS) net- 
works using tetrafunctional crosslinker tetrakisdimethyl- 
siloxysilane. The number- and weight- average molecu- 
lar weights (Mn and Mw, respectively) of the precursor 
PDMS were evaluated to be 46 600 and 89 500, respec- 
tively, by gel permeation chromatography (GPC). In the 
first system, no diluent is added in the whole process; 
the resulting polymer network shall be referred as 100% 
sample. In the second system, 30% of diluent is added in 
the polymerization process and is kept afterwards. The 
resulting system shall be referred as the 70% sample. All 
experiments were done at the temperature of 40^C The 
physical chain densities were estimated to be b.SSmol/m^ 
for the 100% sample and 3.22mol/m^ for the 70% sam- 
ple. [l3| The system prepared by Kawabata is a vulcan- 
ized isoprene rubber (IR) with polydipersed distribution 
of chain length. No structure data is known from the 
preparation method. 



IV. DATA FITTING 

Fitting of stress-strain data using XGR theory is sim- 
ple. From Eq. ()2.2p the stress components is linear 
in /icMi- The variance is therefore a positive-definite 
quadratic function of these two parameters: 



yZ [/io ctq (i) + /ii cri{i) - (Jexp ii)y 



(4.1) 



where the sum runs over all data points. Minimization 
of this variance is a trivial matter, since it is guaranteed 
that there is only one minimum. 

By contrast, fitting of stress-strain data using EV the- 
ory is more complicated. The slip-link theory has four 
parameters, Ci,C2,7?,a, see 12.31 Among these, 77 and 
a enter the theory in an essentially nonlinear way. We 
found that the variance as a function of four parameters 
C'i,C2,?7,a may exhibit multiple local minima that are 
approximately degenerate. In this case, the parameters 
will not be able to be determined using stress-strain data 
alone. For Urayama's systems, however, the parameter 
Ci in Eq. ()2.3p can be determined [2|, [i3| separately using 
structure information. 



A. Urayama's 100% Sample 

We first fit Urayama's data for the 100% sample [2| 
using XGR theory. The best-fitting parameters for /io, Mi 



are very close to each other, as shown in the first row 
of Table [H The fitting of uniaxial stress-strain data is 
shown in Fig. [H The fitting of two stress components 
(Ji and <J2 for general biaxial data are shown in Fig. [2] 
and Fig. [3] respectively. We also perform the data- fitting 
with the constraint /io = Mi • As shown in the second row 
of Table HI the resulting variance is almost the same as 
the unconstrained two parameter fitting. Furthermore, 
the stress strain curves of two cases show no noticeable 
difference. This points to the interesting possibility that 
their theoretical values are indeed the same, for some 
reason to be understood. 
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Let us define an effective chain density n^^^ via 
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ksT, 



(4.2) 



where the temperature T = 40^C = 313K, see 0, ^. 
The best fitting value for /io shown in Table [J then trans- 
lates into an effective chain density n^^^ = 40.0mol/m'^, 
which is much larger than the physical chain density, es- 
timated to be Tie = 5. 887710/ /?Ti^ [2|, [l3]. The parameter 
Mo in XGR theory receives most contributions from en- 
tanglements, rather than from chemical cross-links. 



Theory 


Parameters {MP a) 


Variance (lO^Pa^) 


XGR (mo/ Ml) 


/io = 0.104, /ii =0.112 


2.745 


XGR(/io=Mi) 


Mo = Ml = 0.106 


2.805 


EV (unrestricted) 


(-0.0019,0.165,0.120,0.155) 


2.481 


EV (restricted)^ 


(0.0152,0.150,0.143,0.152) 


2.483 



^The parameter Ci is determined independently. The other three 
parameters are determined by data fitting. 

TABLE I: Best-fitting parameters and variances of Urayama's 100% sample [2|. 
(Ci , C2 , 77, a) respectively. 



The four parameters in EV theory are 



We also fit the same set of data using EV theory. As 
shown inEq. (|T]), an unrestricted fitting with all four pa- 
rameters free to change yields unphysical results: the 
best-fitting parameter mo is negative. Urayama [Tsl inde- 
pendently determined the value of Ci to be 0.0152MPa 
using the knowledge of physical chain density. We have 
also followed the same procedure. The fitting quality thus 
obtained is almost identical to the unrestricted data fit- 
ting. The restricted fitting results of biaxial stress-strain 
data using EV theory are shown in Eq. (|4]) and Eq. (J5]) re- 
spectively. We have checked the individual stress-strain 
data curves, and find no noticeable difference between 
the restricted fitting and the unrestricted fitting. 




FIG. 1: XGR fitting of the uniaxial stress- strain data of 
Urayama's 100% sample 0]. 



B. Urayama's 70% sample 

We have also fitted Urayama's 70% sample using the 
same procedure. This system is prepared in 30% (weight) 
concentration of good solvent. The solvent is kept in 
the measurement stage. As in the 100% sample, the 
best fitting parameters for mo. Mi i^ ^^e XGR theory are 
again very close to each other. The value of mo trans- 
lates into an effective chain density n^^^ ^ 17.5mo//m^, 
which is much larger than the upper bound of real chain 
density ric = 3.22mo//m^, determined independently by 
Urayama [2|, [l3- Hence mo is again dominated by en- 



tanglement effects. The fitting results of uniaxial and 
biaxial stress-strain data are shown in Eq. (|6j), Eq. ([7]) 
and Eq. (J8j) respectively. 

The unrestricted fitting using EV slip-link theory 
yields parameter values that are all physical. The re- 
stricted fitting, with the parameter Ci fixed indepen- 
dently using structure data, yields substantially different 
result for the slippage parameter 77, and a total variance 
that is very close to that of XGR fitting. This is probably 
a generic feature of nonlinear data fitting. 

Caution should be taken in the fitting of either theory 
to stress-strain data of swollen network. Both EV the- 
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FIG. 2: XGR fitting of stress component <7i = 
biaxial deformation of Urayama's 100% sample, 
corresponds to a distinct value of Ai. 
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ory and XGR theory are improvements over the classical 
theory of rubber elasticity, which corresponds to the first 
term of 12.21 and Eq. ()2.3p . The classical theory, how- 
ever, constitutes a zero-th order approximation only for 
un- swollen polymer networks, but not for swollen net- 
works. The first basic assumption, that polymer chains 
obey Gaussian statistics, breaks down when chains are 
swollen by good solvent. It is well known that polymers 
in good solvent obey very different, nonGaussian statis- 
tics [l^. Indeed the fitting quality of both theories for 
the 70% sample becomes noticeably worse than that of 
the 100% sample in the large deformation regime. We 
have tested that the fitting quality becomes even worse 
for samples with higher degree of swelling. 
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FIG. 3: XGR fitting of stress component <72 = df /d\2 for 
biaxial deformation of Urayama's 100% sample. Each symbol 
corresponds to a distinct value of Ai. 



For comparison, we also show the restricted fitting 
results of biaxial stress-strain data using EV theory in 
Eq. dll and Eq. (|5]) respectively. Again no substantial 
difference in fitting quality can be detected between fit- 
tings using two theories. 



Theory 


Parameters {MP a) 


Variance (lO^Pa^) 


XGR(/io//ii) 


/io = 0.045, /ii =0.053 


0.369 


XGR (mo^Mi) 


IJiQ^ IJi^^ 0.048 


0.442 


EV (unrestricted) 


(0.032,0.052,0.419,0.134) 


0.282 


EV (restricted)" 


(0.0083,0.064,0.138,0.155) 


0.354 



•^The parameter Ci is determined independently. The other three 
parameters are determined by data fitting. 



TABLE II: Best-fitting parameters and variances of Urayama's 70% sample [2]. The four parameters in EV theory are 
(Ci, C2, ?7, ol) respectively. The fitting result using EV theory is slightly different from that by Urayama Q]. 



C. Kawabata's Uncharacterized Sample 

We also use XGR and EV to fit Kawabata's data 
I , measured using network with much greater random- 



ness. The results for best-fitting parameters are shown 
in Eq. piip . Again for the two parameter version of XGR 
theory, jjlq^jii are very close to each other. The one pa- 
rameter version of XGR (with /io = /ii) leads to 30% 



a-l (Pa) 



15000O 



100000 




FIG. 4: EV fitting of stress component ai = df/dXi for 
biaxial deformation of Urayama's 100% sample. Each symbol 
corresponds to a distinct value of Ai. 
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FIG. 5: EV fitting of stress component <72 — df/dX2 for 
biaxial deformation of Urayama's 100% sample. Each symbol 
corresponds to a distinct value of Ai. 



the fact that the best-fitting value for the slippage pa- 
rameter ?7 ~ 0.31662 is much larger that of Urayama's 
samples. 

The XGR fit of the uniaxial stress-strain data is shown 
in Eq. (pTj) . There is a slight upturning trend of the data 
in the regime Ai > 2.5, suggesting that the issue of finite 
extensibility of polymer chains is coming to play. This 
is somewhat expected, since the network being analyzed 
has a poly-dispersion of chain lengths. In fact, to avoid 

crl (MPa) 




FIG. 6: XGR fitting of the uniaxial stress-strain data of 
Urayama's 70% sample [2]. 
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more variance, but the stress-strain curves hardly show 
any visual difference. The best fitting parameters for EV 
theory are completely physical. Unlike Uarayama's sam- 
ple, there is no available structure information on Kawa- 
bata's polymer network; the parameter Ci therefore can 
not be determined independently. An unrestricted fit- 
ting using EV theory yields slightly better fitting quality, 
with 25% less variance then XGR theory. The difference 
is however not substantial, given the fact that EV theory 
has two more parameters than XGR. The fitting result of 
Kawabata's data using EV theory is rather robust against 
variations of initial conditions. This is likely related to 



FIG. 7: XGR fitting of stress component ai = df/dXi for 
biaxial deformation of Urayama's 70% sample. Each symbol 
corresponds to a distinct value of Ai. 



dealing with the issue of finite extensibility, we have only 
used biaxial data with Ai < 1.9 for data fitting. [18] The 
biaxial stress-strain data, as well as their fit using XGR, 
are shown in Eq. ([12]) and Eq. (fT3|) . As one can see in 
Eq. ([T2j), the fitting quality becomes systematically worse 
as Ai increases. 
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For comparison, we also show the fitting results of bi- 
axial stress-strain data using EV theory in Eq. ([H]) and 
Eq. ([15]) respectively. Again no visual difference in fit- 



ting quality can be detected between fittings using two 
theories. 



Theory 


Parameters {MP a) 


Variance (lO^Pa^) 


XGR(/io//ii) 


(^o,/ii) = (0.294,0.249) 


2.89 


XGR(/io=/ii) 


/io = Ml = 0.281 


4.07 


EV 


(Ci,C2,ry,a) = (0.22,0.25,0.32,0.12) 


2.13 



TABLE III: Best-fitting parameters and variances of Kawabata sample |3|]. 
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FIG. 8: XGR fitting of stress component a2 = df jdX^ for 
biaxial deformation of Urayama's 100% sample. Each symbol 
corresponds to a distinct value of Ai. 



FIG. 10: Restricted EV fitting of stress component a^ — 
df /d\2 for biaxial deformation of Urayama's 70% sample. 
Each symbol corresponds to a distinct value of Ai. 
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FIG. 9: Restricted EV fitting of stress component ai = 
df/dXi for biaxial deformation of Urayama's 70% sample. 
Each symbol corresponds to a distinct value of Ai. 



FIG. 11: XGR fitting of the uniaxial stress-strain data of 
Kawabata's sample 



V. DISCUSSION AND CONCLUSION 

In this work, we have compared the XGR theory of 
rubber elasticity [l| with the slip-link theory using the 
biaxial stress-strain data of three distinct polymer net- 
works with very different network structures, synthesized 
by Urayama ^ and Kawabata ^] respectively. With two 
adjustable parameters, the XGR theory fits three sets of 
stress-strain data in the whole range of deformation, with 
comparable quality as the slip-link theory, which has four 
tunable parameters. It is probably more remarkable that 
almost the same fitting quality is achieved by the one pa- 



rameter version of XGR. The technical advantages asso- 
ciated with such a one-parameter fitting scheme is rather 
substantial comparing with various entanglement theo- 
ries, which typically contain no less than three parame- 
ters. 

The conceptual differences between XGR and entan- 
glement theories are also worth commenting. Same as the 
classical theory, almost all entanglement theories of rub- 
ber elasticity theory can be categorized as molecular level 
theories, where a many body problem is reduced to that 
of a single polymer chain. By contrast, XGR attributes 
the failure of the classical rubber theory to the entropy 
associated with the collective long wavelength fluctua- 




FIG. 12: XGR fitting of stress component ai = df/dXi 
for biaxial deformation of Kawabata's sample. Each symbol 
corresponds to a distinct value of Ai. 



FIG. 15: EV fitting of stress component a2 = df/dX2 for 
biaxial deformation of Kawabata's sample. Each symbol cor- 
responds to a distinct value of Ai. To improve the visibility 
of the plot, only data with Ai < 1.3 are shown here. 




FIG. 13: XGR fitting of stress component a2 = df/dX2 
for biaxial deformation of Kawabata's sample. Each symbol 
corresponds to a distinct value of Ai . To improve the visibility 
of the plot, only data with Ai < 1.3 are shown here. 
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FIG. 14: EV fitting of stress component ai = df/dXi for 
biaxial deformation of Kawabata's sample. Each symbol cor- 
responds to a distinct value of Ai. 



tions of cross-links, which is manifestly beyond the ca- 
pacity of all molecular level theories. Note, however, the 
strong interplay between long wave length elastic fluctu- 
ations and generalized elasticities has been firmly estab- 
lished in the past years, and actually has been one of the 
main theme of modern condensed matter physics. Fur- 
thermore, XGR demonstrated that the effects of these 
long wavelength fluctuations should also be detectable 
using scattering experiments. By contrast, the effects 
of entanglement on rubbery elasticity in the regime of 
large deformation, though are widely believed and have 
been explored for several decades, have little experimen- 
tal substantiation beyond fittings of stress-strain data. 
In this perspective, it is rather natural to impose the 
following question: what feature in stress-strain data of 
polymer networks is still unexplained after the effects of 
long wave-length elastic fluctuations are taken into ac- 
count? While the answer to this question suggested by 
the present work is " probably nothing" , careful and ex- 
tensive comparison using more sets of data is definitely 
necessary to clarify the situation. 

Perhaps, the most direct experimental test of entangle- 
ment theories would be achieved by a study of polymer 
networks with no entanglement at all. Such a network 
can in principle be prepared by crosslinking a semi-dilute 
polymer solution at c* concentration [l^, where neigh- 
boring polymer coils are barely touching each other, but 
have minimal number of entanglements. After crosslink- 
ing, the solvent should be removed in order to restore the 
incompressibility. Any entanglement theory would pre- 
dict that the elasticity of such a c* network is described 
by the classical theory, while according to XGR theory, 
the elasticity should still be described by the same equa- 
tion Eq. ([22]). 

The current version of XGR theory Eq. (|2.2p only take 
into account the first order perturbative correction of 
elastic fluctuations. On the other hand, the dimension- 



less parameter controlling the perturbative series is given 
hyTa~^fi [1| , which is of order of unity. Therefore higher 
order corrections are expected to be quantitatively im- 
portant, and the good quality of data fitting using first 
order correction Eq. ()2.2p may seem fortuitous. To test 
the validity of XGR theory beyond leading order, we may 
numerically simulate the elasticity of an incompressible 
network in the presence of large elastic fluctuations. This 
shall be done in the near future. 
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questionable is whether entanglements are responsible for 
the failure of the classical theory in the large deformation 
regime. 

[17] We use the result presented in the original paper by Ed- 
wards and Vilgis ^. 

[18] As is well understood, the upturning of the stress-strain 
curve in the large deformation regime is due to the finite 
extensibility of polymer chains. There is no fundamental 
mystery in the understanding of this issue. XGR does not 
take into account the finite extensibility of chains, and 
therefore is not expected to fit the data in the regime 
of large deformations. It is therefore not appropriate to 
compare XGR with the Slip-link theory in this regime. 



